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O ; ABSTRACT 

Numerous numerical studies indicate that dark matter halos show an almost universal radial density 
profile. The origin of the profile is still under debate. We investigate this topic and pay particular 
\ attention to the velocity dispersion profile. To this end we have performed high-resolution simulations 

with two independent codes, ART and Gadget. The radial velocity dispersion can be approximated as 
function of the potential by of = a(<E>/<I > out) K ( < I ) out — <&), where $ ut is the outer potential of the halo. 
For the parameters a and k we find a = 0.29 ± 0.04 and n — 0.41 ± 0.03. We find that the power-law 
■ asymptote a 2 oc $ K is valid out to much larger distances from the halo center than any power asymptote 

for the density profile p oc r~ n . The asymptotic slope n(r — > 0) of the density profile is related to 

OO i 

the exponent k via n = 2k/ (1 + k). Thus the value obtained for k from the available simulation data 
can be used to obtain an estimate of the density profile below presently resolved scales. We predict a 
continuously decreasing n towards the halo center with the asymptotic value n < 0.58 at r = 0. 

en . 

Subject headings: cosmology:theory, dark matter, galaxies: formation, galaxies: structure, methods: 
analytical, methods: numerical 
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1. INTRODUCTION 

The formation of structures by purely gravitationally 
interacting cold dark matter (CDM) is one of the fun- 
damental paradigms in cosmology. The luminous bary- 
onic matter is embedded in dark matter halos. Numerical 
studies of structure formation, which allow only for grav- 
itational interactions, predict the distribution of galaxies 
and clusters of galaxies in excellent agreement with obser- 
vations. Earlier studies had indicated that not only the 
distribution of halos but also their density profiles depend 
on the underlying cosmological model (Quinn et al. 1986; 
Frenk et al. 1988; Dubinski & Carlberg 1991; Crone et al. 
1994). However, subsequent numerical investigations re- 
vealed that the halo density profiles are almost universal 
and depend neither on the mass of the halo (Navarro et al. 
1996, 1997), the initial fluctuation spectrum, nor the un- 
derlying cosmological model. This similarity of dark mat- 
ter halos is complemented by the discovery of a universal 
angular momentum profile (Bullock et al. 2001). 

With the dramatically increased resolution of simula- 
tions the properties of the profiles can now be studied 
over several orders of magnitude in radius. This allows 
a more detailed study of the central slope of the den- 
sity profile. Navarro et al. (1996, 1997) first gave an an- 
alytic approximation for the density profiles. They ob- 
tained p oc r~ n , where the slope depends on the radius 
according to n(r) = (1 + 3r/r s )/(l + r/r s ), r s denotes the 
scale radius of a given halo. This NFW-profilc implies an 
asymptotic slope of n = 1 for the central density profile. 
Other numerical studies produced a significantly steeper 
inner slope, n w 1.5, (Moore et al. 1998, 1999; Ghigna 
et al. 1998, 2000; Fukushige & Makino 2001; Klypin ct al. 
2001). More recently Power et al. (2003) found no asymp- 
totic slope at all but a continuously decreasing slope with 
values down to w 1.2 approaching the innermost radius 
resolved. In order to determine the inner slope of dark 
matter halos in numerical studies more reliably, more ef- 
forts are needed to increase the resolution of simulations. 

In galaxies the motion of the stars is governed by the 
dark matter halo. Measuring rotation curves in the very 
center of a galaxy allows in principle the determination 
of the central profile of the underlying density distribu- 
tion. Low surface brightness galaxies arc almost unbiased 
by any gas and are thus promising candidates for the de- 
termination of the innermost density profile. A 'cuspy', 
n > 1, dark matter core is ruled out by many studies 
(McGaugh & de Blok 1998; de Blok et al. 2001) includ- 
ing models of the inner rotation curve of the Milky Way 
(Binney & Evans 2001). Salucci & Burkert (2000) have 
inferred a constant central density from the kinematics of 
a large sample of spiral galaxies. However, some recent 
investigations of dwarf galaxies suggest that an innermost 



slope of n « 1 is consistent with rotation curve data (van 
den Bosch & Swaters 2001). Lcnsing studies of clusters of 
galaxies indicate that there is a shallow core, with some 
preference of an isothermal profile in the very center (Scitz 
et al. 1998; Gavazzi et al. 2003). Thus, present observa- 
tions favor an inner slope of n « 1 at most and in many 
cases a considerably shallower one. 

Some mechanisms have been proposed which are able to 
erase a central cusp of dark matter halos. El-Zant et al. 
(2001), e.g., suggest that a core may be eliminated by dy- 
namical friction of an initially clumpy gas distribution. 
Alternatively the occurrence of bar instabilities in the cen- 
ter have been discussed (Weinberg & Katz 2002; Sellwood 
2003). Dekel et al. (2003) have argued that a central cusp 
can be avoided by puffing up infalling halos, which could 
occur due to baryonic feedback processes. 

The uncertainty of the central slopes reflects that up to 
now dark matter halo profiles have only been obtained on 
a empirical basis. Efforts are made to calculate the profiles 
analytically and semi-numerical models arew used to cal- 
culate the profile using the growth rate of halos. Since the 
mass accretion depends on the environment those mod- 
els may indicate that the profile should also vary with 
the supposed cosmology (Syer & White 1998; Nusser & 
Sheth 1999; Lokas & Hoffman 2000). An alternative ex- 
planation was given by Taylor & Navarro (2001); Navarro 
(2001). They argued that recurrent merging results in a 
phase-space density profile which decreases according to a 
power-law. As a result they found that the central slope 
should be as low as n — 0.75. One may summarize the 
results about dark matter density profiles by stating the 
open questions. Why do halos show a profile close to the 
NFW-profile, why are the profiles almost universal, and 
what is the central slope? 

In this paper wc investigate the velocity dispersion pro- 
file of dark matter halos that were simulated with high 
resolution and using two independent codes. In Sec. 3 we 
present a suitable formula that approximates the disper- 
sion profile as a function of the potential. The probability 
distribution of the particle velocities within shells of dif- 
ferent radii is investigated in Sec. 4. This allows us to 
interpret the individual factors which are present in the 
suggested formula. The analysis of the possible asymptotic 
behavior of the profiles in Sec. 5 demonstrates that the re- 
lation between velocity dispersion and potential is strongly 
restricted by physical requirements. Moreover, limits for 
the inner slope of the density profile can be set. For consis- 
tency we show in Sec. 6 that the numerical density profiles 
can be recoverd by integrating the basic equations. 

2. RADIAL PROFILES FROM HIGH-RESOLUTION 
SIMULATIONS 
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Fig. 1. — The radial profiles with respect to different quantities 
for two halos (galaxy sized and cluster sized). Upper panel: The 
density p/po (open squares, scaling given by left-hand y-axis), the 
potential (crosses, scaling by right-hand y-axis), the radial ve- 
locity dispersion o> (solid squares, scaling by right-hand y-axis) and 
the anisotropy (3 (triangles, scaling by right-hand y-axis) are plotted. 
Potential and dispersion are given in units of [10 3 km/s] what allows 
us to compare immediately these quantities to the velocity p.d.f. 
Fig. 2. The density is approximated by generalized NFW-profilc 
Eq. (1) (solid line) with parameters r s = 0.42 Mpc, 5 C = 1.01 X 10 4 , 
n = 1.54 and r s = 0.034 Mpc, 8 C = 7.1 X 10 4 , n = 1.11 for the 
cluster-sized halo C16 and the galaxy-sized halo Gal, respectively. 
The anisotropy is approximated by (3 = 0.27 * "I'/'I'out. The values 
of "l^out are taken from Tab. 1. Lower panel: Predictions by the left 
hand side of the Jeans equation (2) including the anisotropy term 
(crosses) and excluding it (open squares) match the circular velocity 
(solid squares). 



We have performed several dark matter high-resolution 
simulations. We have assumed a spatially flat cold dark 
matter model with a cosmological constant favored by 
most current observations. For the first set of simula- 
tions we have used the Adaptive Refinement Tree (ART) 
A^-body code (Kravtsov et al. 1997) and the following cos- 
mological parameters: ACDM with Q,m = 0.3, — 0.7, 
(Ts = 0.9, and h = 0.7. Within a low mass resolution 
simulation (128 3 particles in a 114 Mpc cubed simulation 
box, m part = 2.9 x 10 10 M Q ), we have identified clusters 
of galaxies. From this set we have selected 8 candidates 
with different masses and merging histories and added 5 
smaller clusters/groups for numerical load balance. We 
then re-simulated the clusters with higher mass resolution. 
With particle masses of 4.6 x 10 s M Q a typical cluster and 
its environment contains more than one million particles. 
The highest force resolution with 9 refinement levels was 
0.9 kpc. Subhalos with masses above 4.6 x 10 10 M© are 
well resolved. A typical cluster contains more than 150 
such subhalos. The simulations were done using an MPI 
version of the ART code where each of eight nodes followed 
the evolution of one or two clusters. Another simulation 



within a box of 36 Mpc box length contains a galaxy-sized 
halo. The region containing this halo was simulated with 
an effective resolution of 1024 3 particles, i.e. with a mass 
resolution of 1.7 x 10 6 M©. The highest force resolution 
with 10 refinement levels was 0.1 kpc. 

In addition, we investigate the properties of a cluster- 
sized halo obtained by a high-resolution simulation with 
uniform mass for all particles. This halo is resolved by 
more than one million particles since a huge entire number 
of particles and a comparatively small simulation box was 
used, namely 300 3 particles and a box size of 30 Mpc. The 
initial conditions were set up according to the cosmological 
model: h — 0.65, Om = 0.3, and f^A = 0.7. The simula- 
tion has been performed using the public GADGET-code 
(Springel et al. 2001). For comparison we consider an- 
other simulation, performed with the GADGET-code, with 
the same initial conditions as the cluster-sized halo C16 
simulated with ART. See Tab. 1 for a compilation of the 
halo properties. 

We have determined the radial density profiles for all 
considered objects at redshift z = 0. In Fig. 1 the den- 
sity profiles of a cluster-sized and a galaxy-sized halos are 
shown. Both profiles can be fitted reasonably well by a 
generalized NFW-profile with a free parameter n for the 
inner slope 

p/po = S c (r/r s )- n (l + r/r s y (3 - n \ (1) 
The parameters n, 5 C , and r s are determined for each halo 
by a least-square fit to the mean densities in radial bins 
up to the virial radii. For the halo C16 we obtain an in- 
ner slope of n w 1.54, which corresponds to the Moore- 
profile and which is in agreement with the results given by 
Fukushige & Makino (2001). They analyzed 12 halos with 
various masses and found an inner slope of about n ~ 1.5 
for all of them. In contrast, for most clusters of our sample 
we found smaller inner slopes, see Tab. 1. 

A relaxed spherical halo of collisionless particles is com- 
pletely described by the radial profiles of the density p(r), 
the radial velocity dispersion of = (v r — Tv) 2 , where 
is the mean radial velocity in a spherical shell with mean 
radius r, and the anisotropy of the dispersion (3=1 — 
of/of, where at denotes the tangential velocity disper- 
sion. The potential can be obtained by mass integration 
<& = G J Q r drM(r)/r 2 (Poisson equation), where G de- 
notes the gravitational constant. The radial profiles are 
related by the Jeans equation (Binney & Tremaine 1987) 
d(po1 



dr 



2p 2 d<£ 

— Pcv = -P-J-) 

r ar 



(2) 

which describes a steady-state halo whose particles move 
on collisionless trajectories in a spherical potential self- 
consistently generated by the particle distribution. The 
results of numerical simulations confirm that dark matter 
halos fulfill the Jeans equation at least up to the virial ra- 
dius (Thomas et al. 1998). This leads to the conclusions 
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Fig. 2. — The distribution of the absolute value of the particle ve- 
locities |v| within spherical shells for the galaxy-sized halo Gal. The 
distributions shown in the upper panel are taken in the shell around 
the radius r CTjmax with maximum dispersion and in a second shell 
with smaller radius. The distributions shown in the lower panel are 
taken in shells starting from r CTimax to larger radii. For comparison, 
a Maxwell distribution is shown centered at the maximum of the 
velocity distribution within one of the shells. 

that (i) the halos are relaxed, (ii) the particles in the halos 
are moving on orbits determined by a mean potential, and 
(iii) the small-scale particle-particle interaction is negligi- 
ble, i.e. close encounters are rare. 

We have tested whether our halos fulfill the Jeans equa- 
tion. The circular velocities v c are defined by v\ — GM(r)/r 
r<f>'. Using the Jeans equation they can also be given by 
— r(pa^)'/p — 2/3of , where the prime denotes d/dr. In 
Fig. 1, lower panels, we compare the circular velocities 
obtained by both ways. Furthermore, we compute the 
circular velocities neglecting any anisotropy contribution 
(0 = 0). The differences between the circular velocities at 
given radius are small and, hence, the halos can approx- 
imately be described by the Jeans equation assuming an 
isotropic velocity dispersion. 

3. VELOCITY DISPERSION AS A FUNCTION OF THE 
POTENTIAL 

A spherically symmetric halo with isotropic velocity dis- 
persion can be described using three radial functions, e.g. 
density, potential and velocity dispersion. These profiles 
arc related by the Jeans and the Poisson equation. Thus, 
if one profile or any supplementary relation between a, p 



or $ is given, the two others are determined. An invitingly 
simple way to get the profiles is to assume that the velocity 
dispersion is constant within the virial radius. However, 
the numerical results indicate that the velocity dispersion 
has a pronounced maximum at a radius r = r aimax . For 
example, the galaxy-sized halo Gal shows the maximum 
velocity dispersion of 155km/s at the position r w 30kpc. 
In contrast the dispersion is only 90 km/s at the virial 
radius and to 80 km/s at the innermost (r « lkpc) radius. 

On the basis of the numerical results we make an heuris- 
tic ansatz for the relation between the velocity dispersion 
and the potential at a given radius. Physical arguments 
which lead to this ansatz are discussed in Sec. 4. We ap- 
proximate the dispersion as a function of the potential. 
This can be done because the potential increases monoton- 
ically with radius and, consequently, a one-to-one mapping 
between radius and potential must exist. The data can be 
approximated, see Fig. 3, by the relation 

$ \ K 

3— (*out-*), (3) 

*out / 

where a, n are free parameters and <E> ut is the maximum 
potential reached for large radii. We denote this relation 
between radial velocity dispersion and potential in the fol- 
lowing as VDPR. Since an arbitrary constant value can 
always be added to the potential we have set <&(0) = 0. 
This normalization of the potential is used throughout 
this paper. We perform least-square fits to the data of 
the different halos up to radii of about two times the virial 
radius. Some of the halos have a companion or substruc- 
tures which results in a narrow peak in the dispersion pro- 
file, e.g. halo C110. We exclude those peaks from the 
approximation if they are beyond the virial radius. The 
mean values obtained for the dimensionless parameters are 
a = 0.29 ± 0.04 and n = 0.41 ± 0.03, where the errors are 
standard deviations. 

The primary halo-specific parameter is the outer poten- 
tial $out- Also the parameters a and k vary, but their scat- 
ter is small. It may be caused by substructures still present 
in the halo, imperfect relaxation, or even nearby structures 
which disturb the spherical symmetry. The smallness of 
the scatter may indicate that perfectly relaxed halos have 
only one free parameter, namely the outer potential. Con- 
sequently the VDPR may be considered to be equivalent 
to an equation-of-state since it gives the velocity disper- 
sion as a function of the local potential, independent of 
the radius. On the other hand it incorporates global prop- 
erties, namely the difference of the local potential to the 
central and to the outer potential. This reflects that the 
local particle-particle interaction is small and particles do 
not significantly exchange energy. Therefore, the local ve- 
locity dispersion must depend on the global properties. 
However, if for all perfectly relaxed halos both parameters 
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a and k are the same, the VDPR may serve as an addi- 
tional equation which allows to close the system of Jeans 
and Poisson equation. 

4. VELOCITY DISTRIBUTION IN SPHERICAL SHELLS 

The VDPR consists apparently of two parts. The first 
factor ($/$out) K is dominant in the center of the halo 
whereas the second factor ($ ut — 3>) is dominant in the 
outer region. Also the shape of the velocity probability 
density function (p.d.f.) in a shell with a small radius 
differs significantly from the shape in a shell with a large 
radius. 

In Fig. 2 we provide the velocity p.d.f. determined 
within several spherical shells of the galaxy-sized halo Gal. 
Let us first consider the velocity p.d.f. h(\v\,r) in the up- 
per panel in a shell around 22 < r < 26 kpc close to 
the radius rv imax , where the velocity dispersion reaches its 
maximum. In comparison to a Maxwellian distribution the 
velocity p.d.f. ft(|v|, r ffmax ) is flatter at the small-velocity- 
tail and shows a steep break towards high velocities. The 
second p.d.f. for the shell shown in the upper panel has 
a radius much smaller than rv. max . The velocity p.d.f. is 
much broader compared with a Maxwellian distribution. 
The probability peaks at « lOOkm/s. Albeit there are 
also particles with velocities > 300 km/s - but below the 
escape velocity - the probability to find a particle with 
such a large velocity is very small near the center of the 
halo. This effect could be due to two reasons: first, the an- 
gular momentum of a high velocity particle must be very 
small to get close to the center and second, such a par- 
ticle spends only very little time along its orbit near the 
center. Thus the probability to find particles with a high 
velocity decreases towards the center of the halo, whereas 
the orbit of slow particles is entirely in the vicinity of the 
center. Thus the potential well acts to concentrate slow 
particles in the central regions. Since there is no character- 
istic radius for this concentration process one may expect 
that the velocity dispersion scales at small radii in a self- 
similar manner with the potential, i.e. of oc $ K . This is 
represented by the first factor of the VDPR. 

Let us now consider the shape of the velocity p.d.f. 
within shells around and above the virial radius, see Fig. 2, 
lower panel. With increasing radii of the shells the entire 
distribution shifts to smaller velocities. The shape of the 
velocity p.d.f. itself almost does not change. The proba- 
bility distribution for small velocities, i.e. below the maxi- 
mum of the distribution v max , converges to the Maxwellian 
distribution. For larger velocities the steep break is still 
present. An obvious reason for the shift is that the ki- 
netic energy decreases when a particle moves out of the 
potential well. This implies that the characteristics of the 
velocity p.d.f., e.g. the position of the maximum velocity 




Fig. 3. — The dependence of the radial velocity dispersion / <I>out 
on the potential $/$ ou t. The data obtained from the simulations 
are approximated by the VDPR (3). 

t) max (r) and the velocity dispersion <x(r), depend directly 
on the potential $(r). Particles having velocities different 
from zero for r — > oo, would leave the halo, i.e. they are 
not bound. Consequently the dispersion should vanish for 
large radii for an isolated halo with no additional matter 
around. Thus, interpreting the shift of the velocity p.d.f. 
as a result of the effective potential well and assuming that 
the shape of the p.d.f. is constant for all radii suggests a 
dependence between dispersion and potential according to 
cr 2 (r) cx ($out — < f( r )) at large radii. This is represented 
by the second factor in the VDPR. 

5. ASYMPTOTIC PROFILES 

We now consider the asymptotic behavior of the radial 
profiles at small, r <^r S} and at large radii, r ^> r s . 

First we discuss the outer region, where the dispersion 
can be approximated by of | r3>T . s = a ($ out -$(r)). Using 
the Jeans equation allows we can relate the factor a to 
the slope of density profile. We assume isotropic disper- 
sion, substitute the radial dispersion in the Jeans equation 
by the potential, integrate the Jeans equation and finally 
obtain the relation 

| r>>ra = c ($ out - Syr, 7 = (l-a)/a. (4) 

Since, both, (<f> ut — and the density p, are monoton- 
ically decreasing functions for dark matter halos, the ex- 
ponent must be positive. Therefore, the factor a has to 
be in the range < a < 1. We can constrain a further by 
inserting Eq. (4) into the Poisson equation. This results 
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Code 


^Virial 


-^virial 


-^virial 


c 


n 


$out 


K 


a 






[Mpc] 


[M ] 








[km 2 /s 2 ] 






Gal 


ART 


0.25 


9.2 x 10 11 


1066000 


7.4 


1.11 


1.8 x 10 5 


0.45 


0.33 


CIA 


Gad 


0.90 


7.0 x 10 13 


1031000 


8.8 


1.17 


2.8 x 10 6 


0.44 


0.28 


C16 


ART 


1.52 


1.46 x 10 14 


646000 


3.6 


1.54 


5.2 x 10 6 


0.41 


0.23 


C16b 


Gad 


1.30 


1.34 x 10 14 


590000 


3.6 


1.47 


5.2 x 10 6 


0.41 


0.24 


C17 


ART 


1.27 


1.23 x 10 14 


541000 


5.2 


1.14 


3.8 x 10 6 


0.37 


0.29 


C15 


ART 


1.19 


1.02 x 10 14 


452000 


5.0 


1.12 


3.3 x 10 6 


0.41 


0.33 


C110 


ART 


0.97 


5.59 x 10 13 


247000 


5.7 


1.27 


2.7 x 10 6 


0.38 


0.28 


C19 


ART 


0.90 


4.38 x 10 13 


193000 


8.9 


0.99 


2.4 x 10 6 


0.47 


0.34 


C112 


ART 


0.77 


2.79 x 10 13 


123000 


3.5 


1.46 


1.6 x 10 6 


0.40 


0.29 


Mean 
















0.41 


0.29 



Table 1 

Parameters of the halos obtained by high resolution iV-body simulations. The virial radius is determined according to 
M v ir I (3/4)7rrvi r = 178 fio' 4j p C rit, where p C rit denotes the critical density for a flat universe, see Eke et al. (1998) The density 
profile of the halos is approximated using a least-square fit to the (logarithmic) density values by a generalized NFW-profile 
Eq. (1). The concentration is given by c = r v i T /r s . The relation between the radial velocity dispersion crj: and the gravitational 
potential is approximated by the VDPR (3) using a least-square fits with the free parameters a, k and Bout- 



in the Lane-Emden equation (Binney & Tremaine 1987) 



$) \ + 47rGc($ out - $) 7 = 0. (5) 



One solution of the Lane-Emden equation is a power-law 
for the relative potential 

(*out - *) oc r-"\ m = 2/( 7 - 1) 

and with Eq. (4) also for the density 

Hr»r a OC ($ out - $) 7 OC r~ n , 

In order to obtain a solution for a finite mass halo the 
slope of density profile must be sufficiently steep, namely 
n > 3. Therefore, the factor a must be in the range 

1/4 < a < 1/2. (6) 

Note the asymptotic slope of the NFW-profile, n(r — > 
oo ) = 3, results from a = 1/4. The simulations show that 
the anisotropy of the velocity dispersion increases with ra- 
dius and amounts to (3 « 0.3 at the virial radius. This 
anisotropy affects the allowed parameter range for a: In- 
tegrating the Jeans equation under the condition of a con- 
stant anisotropy parameter (3 leads to 

P(^r)\ r>>rs =c*(^-^rr-^. (7) 

Inserting this relation again into the Poisson equation al- 
lows us to find power-law solutions for the potential and 
the density 

($ 00 -$) C xr- m , m = 2(1 - /3)/( 7 - 1) 



profile is restricted by the demand of a finite halo mass 
even for r —> oo, the parameter a has to be in the given 
narrow range. The halos analyzed here fulfill this condition 
even if they have to be truncated at a radius where the 
ambient, infalling matter becomes dominant. 

Let us now consider the central part of the halos. The 
slope of the inner profile is uncertain from simulations be- 
cause of the lack of force resolution and also because of 
n — 27/(7 — 1) = (2 — 2a)/(l nc K°)? r particle statistics. Therefore, it is still an open 

question whether an inner asymptote for the radial den- 
sity profile exists and whether it is possible to determine 
its power index by numerical simulations. We make the 
weaker assumption that the radial velocity dispersion at 
sufficiently small radii is given by a power-law with respect 
to the potential 



(8) 



In the same manner as described above, we insert the re- 
lation (8) into the Jeans equation, perform the integration 
and obtain 



p($) 



Po 



exp(- 



(fr/iw) 1 - 

a(l — k) 



■). (9) 



where po is a free scaling constant. For k > 1 the assump- 
tion of a power-law with respect to r for $ would result in 
an exponential-like singularity of the density and infinite 
mass at final radii. In order to avoid this, the exponent k 
has to be within the range 

0<k<1. (10) 
^lr>r s x r n = 2(7 ~ ~ 1) — 2(1 — a — a /3) /(I Restricting k in this way the exponential term goes to 

Therefore, the factor a has to be in the range 1/(4 - 2(3) < unity at $ ^ at sufficient small radii, i.e. ($/$ ut) 1_K < 
a < 1/2. Given the positive anisotropy by the simulations, a(l - n) and the inner asymptotic profile for the density 
this narrows the possible parameter range for a. Thus, is given by 

supposing that the outer asymptotic slope of the density p($) | r<r = po ($/$out) _K • (H) 
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Note, for any <I> sufficiently small with respect to <3> ou t 
the asymptote Eq. (8) has the form of a power-law. In 
contrast, the density approaches the asymptotic behavior 
not until ($/$ ut) 1_K "C a(l — k). In particular, if k is 
allowed to get close to unity, this latter condition becomes 
a very strong restriction and is valid only for very little 
potential values. Thus, the power-law like asymptote for 
of can be adopted much further out from the halo center. 
Having this in mind below we will be able to draw con- 
clusions with respect to the density profile near the halo 
center. The power-law ansatz $ oc r m provides a solu- 
tion for the radial profiles. More precisely, by inserting 
the ansatz into the Poisson equation we can determine the 
free constants. For the profiles we obtain 

*(r)l r « r .=*a»t (r/r ) 2/(1+K) 

P(r) \ r<<r = Po (r/ro)" 2K/(1+K) , 



(12) 
(13) 



with 



r 2 - 
r o — 



2ttG (1 + k) 2 po ' (U) 
The free, halo-specific parameters for the inner asymptotes 
are $ out and ro (or po)- In the next section we compute the 
density profiles by integrating the basic equations. The in- 
tegration will be performed starting from the small radius 
ro up to large radii, where the value <fr ut will be reached. 
The inner asymptotes given above allows us to calculate 
the necessary inner boundary conditions at r . 

6. INTEGRATING THE JEANS EQUATION 

Using the VDPR we can compute the profiles for a given 
halo, defined by two parameters, e.g. $ out and r . The 
inner asymptotes are used to set the initial conditions. 
The profiles are obtained by integrating the set of ordinary 
differential equations 



dM 

dr 
M 

dr 

dr 
dp 
dr 
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dr 



d<3> 



+ 



2/3a 2 



-K-lj 



$out, 

dcr 2 

dr ' dr r 
using a Runge-Kutta method. Wc compare the profiles 
obtained in this way with those from the numerical sim- 
ulations. The parameter $ out was already calculated by 
approximating the VDPR, see Tab. 1. We choose ro as 
second free parameter, since it causes only a stretch along 
the radial axis. This parameter can easily be adjusted 
by approximating the maximum position in the dispersion 
profile of(r). First, we assume that the velocity dispersion 
is isotropic, (3 = 0. The resulting density profiles approxi- 
mate roughly the results from the simulations. Clear devi- 
ations occur in the fringes of the halo. This are the regions 




Radius r/r viria i 



Fig. 4. — Integration of the Jeans equation (2). The halos Gal, 
CIA and C16 (triangles, circles and squares, respectively) are in- 
tegrated using the VDPR (3) with the mean values for a and k 
given in Tab. 1. The profiles are integrated for isotropic velocity 
dispersion, j3 = 0, (dashed line) and for an anisotropy according to 
P = 0.27$/$ O ut, see Fig. 1. Initial conditions are calculated from 
Eqs. (12) and (13). The halo-specific values <I>out are taken from 
Tab. 1. For the scale radius ro we obtained from the dispersion 
profile al(r) the values 0.025 Mpc, 0.12 Mpc and 0.13 Mpc for the 
halos Gal, CIA and C16, respectively. 



where considerable anisotropy in the velocity dispersion is 
present. The anisotropy can be reasonably approximated 
by (3 = (So^/^out with Po ~ 0.27. After integration, this 
clearly leads to a better concordance of the density profiles, 
see Fig. 4. For the fit of (3 a dependence on the potential $ 
was assumed in order to avoid the implementation of ad- 
ditional parameters. It can be easily seen that (3 roughly 
follows the shape of $, see Fig. 1. In addition, the results 
are almost insensitive with respect to variations of the de- 
pendence (3 on r as long as [3 is a monotonic function. 

For the halo C16, the integration of the differential equa- 
tions leads still only to a poor approximation of the nu- 
merical results. This halo shows clear signs for undergoing 
still a merging process. The dispersion profile shows a peak 
at w l/8r v ; r , which is caused by a remnant of a merger 
between the cluster and a group about 1 Gyr ago. Thus 
this halo is still in the process of relaxation and therefore, 
can be only roughly approximated by spherical, relaxed 
system. 

The profiles of the halos CIA and Gal are recovered by 
the integration using the mean over all halos for the pa- 
rameters a and k. However, all profiles obtained from 
the simulations are slightly steeper in the center than the 
integrated curves. Dekcl ct al. (2002) found that the in- 
ner profiles of disturbed halos are steepened. Therefore, a 
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more cuspy central profile, compared to a perfectly relaxed 
system, may be caused by some recent merger activity. 

The integrated profiles show a sharp break for radii 
much larger than the virial radius. They do not show the 
asymptotic density profile expected from the discussion in 
Sec. 5. This is due to the fact that the outer potential 
$out is reached at a finite radius, at which the dispersion 
and the density vanishes. However, since the radius where 
the dispersion vanishes is much larger than the virial ra- 
dius of the considered halos, this outer slope cannot be 
determined by simulations designed to investigate struc- 
ture formation. 

7. DISCUSSION AND CONCLUSION 

The numerous numerical studies of the structure of dark 
matter halos indicate that perfectly virialized systems show 
presumably a universal density profile. The latter can be 
approximated reasonably well by the NFW-profile or a 
generalized version of it. Similarly we have introduced the 
relation between velocity dispersion and potential Eq. (3) 
(VDPR): it characterizes the profile by very few param- 
eters. In contrast to the approximations of the density 
profiles the VDPR is implicit in the sense that it as a 
function of the potential instead of the radius. However, 
by introducing the VDPR we can put restrictions even on 
the inner asymptotic slope of, both, the velocity dispersion 
profile and the density profile. The data can be approxi- 
mated very well by the newly introduced VDPR. Almost 
all halos clearly show a power-law for of (<&) cx $ K near 
the center, the best resolved halos do so over one order of 
magnitude. We suppose that the obtained exponent k re- 
flects also the inner asymptotic slope of of (<j>). As a result, 
the innermost density profile power index n is related to 
k by n — 2k/ (1 + k). Using the mean value k — 0.41, de- 
termined from all presented simulations we get n = 0.58. 
The obtained standard deviation An = 0.03 permits only 
a negligible variation of the exponent n. Thus, our pre- 
diction for n is below the theoretical prediction by Taylor 
& Navarro (2001). They argued that the inner asymp- 
totic slope is given by n — 0.75. Our result also indicates 
that the innermost slope of the density profiles is less steep 
than that obtained from the density profiles directly. This 
is not a discrepancy since the prediction n = 0.58 is valid 
for scales smaller than those which are currently resolved 
by numerical simulations. Eq. (9) indicates that the den- 
sity profile may steepened by an exponential term in a 
range where the velocity dispersion profile already follows 
a power-law. The density profile only becomes a power- 
law if ($/$ out ) 1_K <C o(l — k) is fulfilled. For instance, the 
impact of the exponential term becomes smaller than 10% 
if the potential is below 0.001 x $ out , using the obtained 
values for n and a. 



In a theoretical analysis of all possible asymptotes Miicket 
& Hocft (2003) found that n should be in the range < 
k < 1/3. Therefore, the inner asymptotic slope of the den- 
sity profile should be as low as n < 0.5. Our numerical 
results are close to this theoretical predictions. Possible 
reasons for the 20% higher value of k could be too strong 
assumptions in the theoretical analysis or steepening of 
the central density profile during relaxation. In fact, it 
is still uncertain to which degree dark matter halos are 
actually relaxed. During the cosmological evolution halos 
merge frequently and show almost always substructures. 
Therefore, the present substructure may have impact on 
the resulting profiles. Nevertheless it is assumed that viri- 
alization has occured to a high degree and the halo con- 
figuration obeys the Jeans equation. The sample of halos 
we considered contains very relaxed halos, e.g. the halo 
Gal, as well as halos with clear deviations from relaxation, 
e.g. the halo C16. Dekel et al. (2002) found that ongo- 
ing merging enforces a central cusp, consequently we may 
have overestimated the mean value of k to some extent. 

Our numerical results do not indicate any significant 
variation of the exponent k towards the center, see Fig. 3. 
Only in the case of the highest resolved halo Gal a trend 
of profile flattening can be noticed. Generally, a core 
with constant velocity dispersion might exist one a sub- 
resolution scale. The existence of such a core would weaken 
the restrictions for k. However, the apparent inner con- 
stancy of k in the numerical simulations rather points to 
the existence of a power-law asymptote. 

A further benefit of analyzing the VDPR is that this re- 
lation enables the determination of the leading processes 
responsible for the quite universal halo density profiles. 
The VDPR consists of two parts, the power-law $ K and 
the relative potential $ ut — The latter can be under- 
stood as the gain in energy of the particles falling into 
the potential well. Utilizing the virial theorem 2T(r) = 
U(r), where T(r) = (1/2) £ r<7 , <r+Ar m iV f and U(r) = 
J2r<n<r+Ar m i(®out ~ $) and assuming isotropy for the 
velocity dispersion, 3of = vf we obtain of = (l/3)($ ou t — 
$). The factor 1/3 is very close to the value a = 0.29. 
Dark matter halos, therefore, can be considered to be lo- 
cally virialized in their outer region. This leads to the 
observed, steep outer profiles. In the inner region of the 
halo the additional factor $ K occurs, which may be at- 
tributed to the motion of relatively slow particles concen- 
trated within the central potential well. Consequently, 
slow particles are more likely to be found in the center. 
This argument is scale-independent and should therefore 
result either in a power-law for of ($) or in a dispersion 
core in the center of the halo. The characteristic and uni- 
versal density profiles are a result of these two mechanisms, 
which do neither depend on the history of the halo forma- 
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tion nor on the cosmological environment. For a perfectly 
relaxed system an universal density profile should be ex- 
pected. 

In summary, we have shown that the radial velocity 
dispersion can be approximated well by the VDPR as a 
function of the potential. We have provided physical argu- 
ments for the analytical form of the VDPR. The numerical 
data already obey the inner power-law over one order in 
magnitude. Considering the inner asymptotics, we have 
shown that the density profile becomes a power-law if the 
potential is below « 0.001 x $ ut- We conclude that the 
density profile can reach its asymptotic slope n < 0.58 
only at scales smaller than those presently resolved. 
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